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ABSTRACT 

We study the absorption lines present in the spectra of the long-period cataclysmic 
variable GK Per during its quiescent state, which are associated with the secondary 
star. By comparing quiescent data with outburst spectra we infer that the donor star 
appears identical during the two states and the inner face of the secondary star is 
not noticeably irradiated by flux from the accreting regions. We obtain new values for 
the radial velocity semi-amplitude of the secondary star, Kk = 120.5 ± 0.7 kms -1 , 
a projected rotational velocity, VkSuii = 61.5 ± 11.8 kms -1 and consequently a 
measurement of the stellar mass ratio of GK Per, q = Mk/Mwd = 0.55 ± 0.21. The 
inferred white dwarf radial velocities are greater than those measured traditionally 
using the wings of Doppler-broadened emission lines suspected to originate in an 
accretion disk, highlighting the unsuitability of emission lines for mass determinations 
in cataclysmic variables. We determine mass limits for both components in the binary, 
M K > 0.48 ± 0.32M Q and M WD > 0.87 ± O.24M . 

Key words: 

binaries: close - stars: cataclysmic variables - stars: fundamental parameters - stars: 
individual: GK Per. 



1 INTRODUCTION 

GK Per (Nova Per 1901; Campbell 1903), is a cataclysmic 
variable (CV) that undergoes dwarf nova outbursts every 2- 
3 years, brightening from 13th to 10th magnitude (Sabbadin 
& Bianchini 1983). Its spectrum is a composite of broad 
emission lines and narrow absorption lines that move sinu- 
soidally in wavelength over the orbital cycle, but 180 degrees 
out of phase, indicating origins associated with both stars. 
The absorption lines are generally photospheric and from 
the secondary star (Kraft, 1964), while the emission lines 
are probably from an accretion disc around the compact ob- 
ject (we explore this further in a future paper). Within the 
sample of known CVs, GK Per has the longest orbital pe- 
riod, P orb = 1.996803 d(0.000007), (Crampton, Cowley & 
Fisher 1986). This makes observing the orbital cycle from 
the ground difficult unless observations are separated over 
several observing seasons or a number of observatories are 
employed contemporaneously at varying longitudes. 

Previously, (Morales-Rueda, Still & Roche 1996, 1999; 
hereafter MSR96 and MSR99), we have presented spec- 
trophotometry observations of GK Per taken during its 
1996 outburst (Mattei et al.1996). In the current paper we 



present spectroscopic observations of GK Per during its qui- 
escence state and focus mainly on the behaviour of the sec- 
ondary star absorption lines, searching for signs of heating 
on the inner surface of the companion. Furthermore we mea- 
sure the radial velocity of the secondary star and combine 
this with previously published velocities. With all the data, 
covering from 1964 to 1995, we calculate a new ephemeris 
for GK Per and a more accurate value for the radial veloc- 
ity semi-amplitude of the companion star, Kk. We measure 
Vk sin i and consequently the mass ratio q for GK Per. A 
lower limit for the white dwarf mass is derived from these 
values which is significantly different to the mass inferred by 
emission line measurements. This provides a further exam- 
ple of the general unsuitability of emission line for measuring 
stellar masses in CVs (Dhillon, Marsh & Jones, 1991). 



2 OBSERVATIONS 

Spectra of GK Per during its quiescent state were acquired 
in 1995, from October 11 to 13, and from October 31 to 
November 2. The first set of spectra was taken at the 2.5 m 



© 0000 RAS 



2 L. Morales- Rueda et al. 



Table 1. Journal of observations. The orbital phases were calcu- 
lated using the ephemeris obtained in Section 3.2. Phase cor- 
responds to superior conjunction of the white dwarf. Labels (1) 
and (2) indicate which observation campaign the data were taken 
at: (1) INT, (2) McDonald. 



Date 


Start 


End 


Start 


End 


No. of 


(1995) 


CUT) 


(Orbital phase) 


spectra 


11 Oct (1) 


22.27 


4.58 


0.027 


0.158 


129 


12 Oct (1) 


22.91 


6.48 


0.541 


0.699 


165 


13 Oct (1) 


3.62 


6.59 


0.140 


0.202 


57 


31 Oct (2) 


5.76 


9.67 


0.199 


0.281 


30 


1 Nov (2) 


5.22 


11.71 


0.689 


0.824 


51 


2 Nov (2) 


4.88 


8.96 


0.183 


0.268 


32 



Isaac Newton Telescope (INT) on La Palma using the In- 
termediate Dispersion Spectrograph. The standard readout 
mode was used with the Tektronix CCD TEK4 windowed to 
1024 x 340 pixels. Exposure times were 35 s with 90 s dead 
time, and the spectral resolution at H/3 was 110 kms -1 . 
The R1200B grating with 1200 lines mm -1 covered the wave- 
length range from AA4136 - A5000 A. 351 spectra of GK Per 
and 8 spectra of K-type template stars were obtained during 
this observing run. Several bias and Hatfield frames were ob- 
tained each night to correct for the bias level and the pixel- 
to-pixel response variations of the chip, and a CuAr frame 
was taken every 15 frames of the target to calibrate the 
spectra in wavelength. A wide slit exposure of the standard 
HD19445 (Oke & Gunn, 1983) was taken every night to cali- 
brate the detector's larger scale response variations. A spec- 
trograph slit orientation of PA 249.1° allowed a 15th magni- 
tude nearby star approximately 25 arcsec ENE of GK Per to 
be employed as calibration for light losses on the slit. A wide 
slit spectrum of GK Per and the comparison star was also 
taken to calibrate the absolute flux scale. After debiasing 
and flat-fielding the frames, spectral extraction proceeded 
according to the optimal algorithm of Home (1986). 

The second set of spectra was taken using the low-to- 
moderate resolution spectrograph ES2 mounted on the 2.1 m 
telescope at McDonald Observatory in Texas. The Til CCD 
camera and a grating ruled at 1200 lines mm -1 were used, 
covering wavelengths AA4196 - A4894 A. Exposure times 
were 350 s with 40 s readout time, and the resolution at 
H/3 was 130 kms -1 . Bias and Hatfield frames were taken 
each night with regularly-spaced arc lamp exposures. The 
standard Feige 110 (Stone, 1977) was observed to flux cal- 
ibrate the data. Simultaneous photometry of GK Per was 
taken during the campaign at the McDonald 36-inch tele- 
scope with the Stiening Photometer (Zhang et al., 1991) us- 
ing UBVR filters (see Skidmore et al. 1997, for a description 
of the filters). The exposure times for the photometry were 

1 s. The combination of the spectroscopy and the photome- 
try yielded 113 spectra of GK Per calibrated in flux. Spectra 
of 7 K-type template stars were also obtained. Tables 1 and 

2 give a journal of the GK Per observations and a list of 
the K-type templates observed during both campaigns and 
their spectral classes respectively. 



3 RESULTS 



Table 2. K star templates observed, with their corresponding 
spectral classification, (a) Houk & Smith-Moore (1988); (b) Nas- 
sau & van Albada (1947); (c) Yoss (1961); (d) Roman (1952); 
(e) Georgelin (1967), (f) Griffin & Redman (I960), (g) Hoffleit & 
Warren (1991). 

Name SpType Ref. Name SpType Ref. 

October 1995 November 1995 



HR190 


Kim 


a 


13 Lac 


KOin 


d. 


f 


HR8688 


Kim 


b 


1 Peg 


Kim 


d. 


f 


HR8415 


K2m 


c 


69 Aql 


K2m 


d 




HR8632 


K3m 


d 


39 Cyg 


K3m 


d 




HR8974 


Kliv 


d 


33 Vul 


K3.5m 


g 




HR8881 


Klv 


c 


3r;Cep 


KOiv 


d. 


f 


HR222 


K2v 


f 


HD197964 


Kliv 


g 




HR8832 


K3v 


d, f 











3.1 Spectral classification of the donor star 

In 1964 Kraft noticed the presence of absorption features 
characteristic of K-type stars in the optical spectrum of 
GK Per. Kraft (1964), Gallagher & Oinas (1974), Cramp- 
ton et al. (1986) and Reinsch (1994) carried out temporal 
studies of these absorption lines and concluded that their 
origin was the mass-donating secondary star of the binary 
system. These authors classified the donor star by compar- 
ing its absorption features with the same features seen in 
K-type template stars. Kraft (1964), Gallagher & Oinas 
(1974) and Reinsch (1994) obtained consistent classifica- 
tions; K2ivp, K2ivp - K2v, and K2v or K3v respectively. 
Kraft (1964) and Gallagher & Oinas (1974) noted changes 
in the line ratios from spectrum to spectrum throughout 
their observing campaign and the classification they gave 
was the mean spectral type of all the observations. The spec- 
tral type of the secondary as measured by Crampton et al. 
(1986) (KOin - KOiv) is not consistent with those of the 
authors mentioned above. Moreover, Crampton et al.(1986) 
conclude that phase-related spectral changes seen by previ- 
ous authors are not present in their own data. 

The spectra of GK Per and the K-type templates pre- 
sented here were binned into constant velocity bins. We used 
the fit to the orbital radial velocity of the secondary calcu- 
lated in Section 3.2 (systemic velocity 7=40.8 kms -1 , ra- 
dial velocity semi-amplitude of the companion Kk = 120.5 
kms -1 ) to shift out the orbital motion of the absorption 
lines with a quadratic rebinning algorithm. The GK Per 
spectra from each night were then binned into two orbital 
phase intervals to investigate the possible changes in spec- 
tral and luminosity classes of the secondary at different 
orbital phases. We selected the regions where there were 
mostly-uncontaminated absorption features, i.e. regions con- 
taining obvious emission lines (AA4265-4399A, AA4435- 
4484A, AA4503-4553A, AA4620-4727A, AA4842-4878A, and 
AA4910-4935A) were masked out from the data. 

In order to classify the secondary star, an optimal sub- 
traction method (Marsh, Robinson & Wood 1994) was ap- 
plied. The K star templates are individually subtracted from 
the total spectrum. The residual is compared to a smoothed 
version of itself, produced by convolving the residual with a 
Gaussian of FWHM = 13 A using the \ 2 statistic. For each K 
star, this procedure is iterative, where the template is scaled 
by a multiplicative free parameter which is determined by 
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Figure 1. The top spectrum is an average of the spectra of GK Per obtained during the first night of observations on the October 
1995 campaign at the INT. The middle spectrum corresponds to the Kliv template HD197964 multiplied by a constant, 0.3xl0 — 3 . 
The bottom spectrum is the residual resulting from the subtraction of the template from the averaged GK Per spectrum and probably 
resembles the spectrum of the accretion flow during quiescence. Some absorption and emission lines present in the spectra have been 
labelled. Most of the absorption lines were identified by Crampton et al. (1986). 



minimising \ 2 - The assumption is that the correctly-scaled 
template spectrum is a good representation of the secondary 
star and the most intrinsically-smooth residual is a good 
representation of the accretion flow spectrum. The tem- 
plate providing the smallest minimum-x 2 is considered the 
best match to the secondary star spectrum. A plot of the 
minimum-^ 2 per degree of freedom found for each template 
for the two orbital phase intervals observed each night can 
be found in Fig. ^| 

The best-fit template is a KllV star (HD197964) in all 
cases. The companion star contributes 69 percent of the total 
light at 4550A, in contrast with 13 percent measured dur- 
ing outburst (MSR99). Fig. 1 presents the total spectrum 
of GK Per before (top panel) and after (bottom panel) the 
subtraction of the scaled Kliv star template (middle panel). 
Identified emission and absorption lines are labelled. This 
spectral classification of the secondary is close to previous 
studies of GK Per during quiescence. This result does not fit 
within the spectral class-orbital period relation fit by Smith 
& Dhillon (1998). This may indicate that the secondary star 
in GK Per is not a main-sequence star, which is not surpris- 
ing as basic evolutionary theories tell us that one should not 



find a main-sequence star in a CV with an orbital period 
much longer than ~9 hours (Patterson 1984). 

The INT data covers approximately the same orbital 
phases as the 1996 outburst data analysed in MSR96 and 
MSR99. The template star nearest in spectral class to the 
donor during outburst at phases, cj> ~ 0.14-0.18 and <j) ~ 
0.64-0.68 is also a Kliv star (MSR99). The secondary star 
shows the same spectral class during both outburst and qui- 
escence at those particular phases. Furthermore, there is no 
convincing evidence for the rapid spectral changes reported 
by Kraft (1964) and Gallagher & Oinas (1974). 



3.2 Radial velocity of the secondary star 

The continuum was subtracted from individual spectra by 
masking wavelength regions with major emission lines, fit- 
ting a spline with 3 knots and subtracting the fits from 
the GK Per spectra. For the template, a two-knot spline 
fit was subtracted. The spectra were binned in velocity to 
ensure that all covered the same wavelength region and had 
the same logarithmic dispersion. Emission line regions re- 
mained masked out and the absorption lines of GK Per were 
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cross-correlated (Tonry & Davis 1979) against those of the 
Kliv template. The template spectrum was broadened by 
the value of Vk sin i obtained in Section. 3.4 before cross- 
correlation was carried out. This is an iterative process be- 
cause to measure Vk sin i we need to know Kk . We carried 
out three iterations after which the values of Kk and Vk sin 
i were stable. 

Cross-correlation of the absorption lines of GK Per and 
the broadened template yielded the radial velocities of the 
secondary star in the binary, relative to that of the Kliv 
template used. The radial velocities obtained were then 
shifted to account for the radial velocity of HD197964; 
-6.5 kms" 1 (Evans 1979). 

We combined these results with previously published ra- 
dial velocities (Kraft 1964, Crampton et al. 1986, and Rein- 
sch 1994), all obtained during quiescence and carried out 
a search for periodicities using the Lomb-Scargle algorithm 
(Scargle 1982). Based on this period search we obtained a 
new ephemeris for GK Per: 

Tq(HJD) = 2450022.3465(0.0006) + 1.9968(0.0008)£ (1) 

We plot the radial velocities folded using this ephemeris 
in Fig. We did not include the radial velocities measured 
during outburst because they show increases and decreases 



steeper than those seen during quiescence (see MSR99 for a 
comparison between radial velocities measured during qui- 
escence and outburst). Although we do not have an explana- 
tion for this behaviour, we decided not to include the out- 
burst data for consistency. The errors on individual mea- 
surements previous to this study were assumed to be equal 
to the mean error of 20 kms" 1 (Reinsch 1994). We fitted all 
the combined data with a sinusoid: 



V — 7 + Kk sin 2tt[4> — 4>o] 



(2) 



and plot it with a solid line in Fig. The best fit pro- 
vides 7= 37.9 ± 0.3 kms" 1 , and K K = 119.3 ± 0.6 kms" 1 . 
For comparison, we fitted only the radial velocities previ- 
ous to this work with a circular function obtaining 7= 22.5 
± 1.7, and K K = 129.0 ± 2.7 kms" 1 . The uncertainties of 
the fit parameters from the combined data are smaller than 
those of the published data, and the two solutions for Kk are 
consistent within 3-er. We also fitted a sinusoidal to our data 
only and obtained 7= 40.8±0.7 kms" 1 and K K = 120.5 ± 
0.7 kms" 1 . This value for Kk is consistent within 1-cr with 
that obtained from the combined data. For clarity we plot 
two orbits in Fig. |^. Our data is plotted only once so it is 
easier to see how good the fit is for previously published 
data. 
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Figure 3. Radial velocity curve obtained by cross-correlating the absorption features of the companion star for the quiescent data sets 
against suitable templates. The radial velocities measured by Kraft (1964), Crampton et al.(1986), and Reinsch (1994) are also plotted. 
The solid line is the circular fit to all the data. The dashed line is the fit to the data previous to this study. The dotted line is a fit to 
the data presented in this paper. Note that although two orbits are plotted, our data is only plotted once for clarity. 



We noticed that between orbital phases 0.4 - 0.8 the 
previously published data shows large deviations from the fit 
to all the data. This is not surprising considering that the ra- 
dial velocity measurements were done using different radial 
velocity standards, of different spectral types in some cases. 
We mentioned before that Kraft (1964) had noticed spec- 
tral changes in the secondary during his observations which 
could have been caused by irradiation of the secondary. This 
would explain the deviations from the fit observed. For this 
reason we have decided to adopt as the best solution the 
fit to only our data, as it was obtained from a consistent 
dataset. 

A second-order sinusoidal fit has been found in the past 
to provide a statistically more-significant fit than a circular 
one. This results from the irradiation of a fraction of the 
stellar surface (Friend et al. 1990). Such a fit in this case 
provides no significant eccentricity, verifying we have not 
detected an irradiated atmosphere in this data. 



3.3 Emission line velocities 

Crampton et al. (1986) measured the radial velocity semi- 
amplitudes of the wings of H/3 to be 34 it 5 kms -1 and 



combined with their absorption line measurements reported 
a mass ratio q (Mk/Mwd) of 0.28, where JV/wd and Mk are 
the masses of the primary and secondary star respectively. 
Although they used only H/3 because the variations in the 
radial velocities for the other emission lines showed more 
scatter, it is clear from their Fig. 3 that even the scatter in 
H/3 gives a very unreliable fit. Adopting the emission line ve- 
locity to be the same as the white dwarf requires the assump- 
tion that the disc emits axisymmetrically, however there are 
many examples where this is not the case (Kaitchuck et al. 
1994). Axisymmetry can generally be minimised by sam- 
pling the velocities of the high velocity line wings only since 
the inner disk is expected to be more axisymmetric than the 
outer (Marsh 1988) and we apply this method to GK Per. 

We use the double-gaussian convolution technique of 
Schneider & Young (1980). This convolves each continuum- 
subtracted, logarithmically-rebinned spectrum with two 
Gaussians of identical width, FWHM = 200 kms" 1 , and 
of separations varying systematically between 400 and 
3000 kms" 1 . The resulting radial velocities were fitted us- 
ing: 

V = 7 - K cm sin 27r[0 - <f> ] (3) 
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Figure 5. The radial velocity curves for H/3 and H7 from mea- 
surements with a double Gaussian separation of 900 kms — 1 . 
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Figure 4. Diagnostic diagram of the Balmer lines, showing 
the radial velocity semiamplitude K em and its fractional error, 
<rx/K cm , as well as the systemic velocity, 7, and the zero phase, 
4>o- The zero phase is relative to the ephemeris derived in Sec- 
tion 3.2. The dashed and dash-dotted lines represent the results 
obtained after applying the Schneider & Young (1980) method to 
H/3 and H7 respectively whereas the solid line presents the results 
for the sum of both lines. 



and the parameters of the fit are plotted in the diagnos- 
tic diagram (Shatter, Szkody & Thorstensen 1986) shown 
in Fig. Faster gas velocities are sampled as the Gaus- 
sian separation increases. The value of K em depends sensi- 
tively on the segment of line that is sampled. Commonly, 
Kwd is considered to coincide with K em when K cm is stable 
over a range of Gaussian separations, i.e. 36 ± 13 kms -1 
in this case, or at the point just before o\R-/K cm increases 
sharply, i.e. at a Gaussian separation of ~900 kms" 1 wich 
corresponds to ~40 kms -1 in this instance. The actual ra- 
dial velocities and the fits obtained are plotted in Fig. ^ 
for measurements taken when the Gaussian separation is 
900 kms -1 . Another approach to the diagnostic diagram is 
the light centres method (Marsh 1988) in which we plot the 
radial velocities measured in velocity space, extrapolate the 
line of points obtained to the K y axis and read off the value. 
Using this method we find K cm ~ 30 kms -1 for H/3 and 
kms -1 for H7 (see Fig. |). 

However, both approaches are not particularly justified 



H/3 



(km/s) 



Figure 6. The light centres of the disc deduced from the radial 
velocity fits of H/3 and H7. The position moves from left to right 
as the Gaussian separation increases which is equivalent to mea- 
suring the velocities closer to the white dwarf. For H/3 the Gaus- 
sian separations go from 400 to 1400 kms -1 in 100 kms' 1 steps 
whereas for H7 they go from 600 to 1200 kms -1 in 100 kms 
steps. In H/3 we have ignored separations above 900 kms -1 in the 
extrapolation to the K a axis as they start departing from K^=0. 



and we find it impossible to determine Kwd from these mea- 
surements. 



3.4 Rotational line broadening 

Fortunately, if we make the assumption that the secondary 
stars rotation is tidally-locked to the orbital period we can 
determine the mass ratio of the binary using only the pro- 
jected radial velocity semi-amplitude and rotational velocity 
of the secondary (Section B.5). This rotational quantity can 



be measured from the width of the absorption lines. Calcula- 
tion of the masses of the components of binary systems from 
Vk sin i measurements is considered a more accurate method 
than measuring the emission line velocities, although Bleach 
et al. (2000) provide some caveats. 

We used the radial velocity fit obtained in Section 3.2 
to subtract the orbital velocity of the absorption features 
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in GK Per, and averaged the continuum subtracted spectra 
for each night. We convolved the Kliv template spectrum 
with a broadening function in which we include a linear limb- 
darkening coefficient of 0.53 (Gray 1992; van Hamme, 1993), 
obtaining a series of templates corresponding to Vk sin i = 
5-150 kms" 1 . Then we applied the same optimal subtrac- 
tion technique as in Section 3.1 to find which broadened 
templates best matched the lines in the average spectrum. 
This technique has been applied by many authors to close bi- 
nary systems (e.g. Southwell et al., 1995; Marsh et al., 1994; 
Drew, Jones & Woods, 1993; Friend et al., 1990). 

The INT and McDonald spectra yielded V K sin i = 64.3 
± 11.0 kms -1 , and 58.7 ± 14.2 kms -1 respectively. We 
adopt Vk sini = 61.5 ± 11.8 kms -1 , the mean of the values 
found for the quiescent data. 



dwarf mass Mwd = 1.44M©. The grey-shaded regions give 
the possible solutions to the mass of the system (see figure 
caption for more details). 



4 CONCLUSION 

We have measured the mass ratio of GK Per, q = 0.55 ±0.21 
using the width and orbital velocity of absorption lines from 
the secondary star, and determine lower limits for the stel- 
lar masses, M W d > 0.87 ±0.24 and M K > 0.48 ± 0.32 M G , 
without using dubious measurements of Kwd ■ The spectral 
type of the secondary star remains constant over the orbital 
cycle indicating that the inner face of the companion is not 
significantly irradiated by accretion radiation. 



3.5 Stellar masses 
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From Kk= 120.5 ± 0.7 kms , measured in Section 3.2, we 
can obtain a value for the mass function of white dwarf of 
0.362 ± 0.006 M Q using: 



Mwd sin 3 i 
(1±Q) 2 



KkP 
2-kG 



(4) 



In Section 3.5 we defined the stellar mass ratio q, where 
M K .Kwd 



M W d 



(5) 



The rotational, Vk sin i, and radial velocity semi-amplitude, 
Kk, of the secondary star are related by 



Vk sini = K K (l + q)R,L K (q) 



(6) 



where Ri JK (q) is the spherically-averaged Roche lobe ra- 
dius of the companion star, scaled to the stellar separation. 
Eggleton (1983) has supplied an analytical approximation 
for this, correct to ~ 1%: 



Rh K (q) 



0A9q 



2/3 



O.i 



< q < OO 



(7) 



7 2 / 3 ± ln(l ± q 1 / 3 ) 

Consequently q = 0.55 ±0.21. This value does not agree with 
q — 0.28 measured by combining absorption and emission 
line velocities by Crampton et al. (1986) and the implied 
velocity semi-amplitude for the white dwarf is Kwd= 66.3 
± 25.3 kms -1 . 

Substituting q into Eq. 4, we obtain a lower limit for 
the mass of the white dwarf of M W d = 0.87 ± O.24M . 
Combining this value with q we then have a lower limit on 
M K of 0.48 ± 0.32 M Q . 

Reinsch (1994) observed GK Per at orbital phases dur- 
ing which potential eclipses would happen but did not mea- 
sure any eclipse effects in the line wings. This implies that 
the white dwarf and inner accretion region are not eclipsed. 
Eclipses of the compact object and inner accretion regions 
require an orbital inclination i > cos -1 Rl k , where Rl k is 
the spherically-averaged Roche lobe radius of the secondary 
star, scaled to the stellar separation, defined in Eq. 6. 

Assuming a value for q of 0.55 ± 0.21, we obtain that 
eclipses would occur if i was greater than 73° . Fig. 7 repre- 
sents the possible values for the masses of the components of 
the system depending on the mass ratio q and the inclination 
i of the binary. Solutions are presented between the mass 
function value of 0.362Mq , and the Chandrasekhar white 
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Figure 7. Possible values for Mk and Mwd from the calculated mass function. Since the system does not show eclipses, i must be 
smaller than 73° (limit obtained from the maximum value of q we calculate). The minimum value possible for Mk found from this work is 
0.48 ± 0.32 Mq. This implies a minimum value for Mwd of 0.87±0.24 Mq. The maximum possible value for Mwd is the Chandrasekhar 
limit (1.44 Mq). The light grey-shaded region gives the possible mass solutions for GK Per. The limits of this region are the minimum 
and maximum possible value of Mwd and the minimum and maximum values of q found from this work . The dark grey-shaded region 
includes the limit imposed by the inclination of the system. Several mass ratio curves are also plotted. 
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